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ABSTRACT 

We present a phase-connected timing solution for the nearby isolated neutron star RX J1308. 6+2127 
(RBS 1223). From dedicated Chandra observations as well as archival Chandra and XMM-Newton 
data spanning a period of five years, we demonstrate that the 10.31-sec pulsations are slowing down 
steadily at a rate of P = 1.120(3) x 1CU 13 s s -1 . Under the assumption that this is due to magnetic 
dipole torques, we infer a characteristic age of 1.5 Myr and a magnetic field strength of 3.4 x 10 13 G. As 
with RX J0720.4— 3125, the only other radio-quiet thermally emitting isolated neutron star for which 
a timing solution has been derived, the field strength is roughly consistent with what was inferred 
earlier from the presence of a strong absorption feature in its X-ray spectrum. Furthermore, for both 
sources the characteristic age is in excess of the cooling age inferred from standard cooling models. 
The sources differ, however, in their timing noise: while RX J0720.4— 3125 showed considerable timing 
noise, RX J1308. 6+2127 appears relatively stable. 

Subject headings: pulsars: individual (RX J1308. 6+2127) — stars: neutron — X-rays: stars 
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1. INTRODUCTION 

The sample of nearby, radio-quiet, isolated neutron 
stars discovered by ROSAT (for a review, e.g., lHaberll 
l2004j) is of particular interest because of the unambigu- 
ous presence of strong, broad absorption features in the 
X-ray spectra of most sources. Since the spectra ap- 
pear thermal, these features almost certainly arise in 
the neutron-star atmospheres. The features have usu- 
ally been interpreted under the assumption of a pure hy- 
drogen composition (not unlikely, given the short set- 
tling times for neutron stars) and a strong magnetic 
field (as suggested by the long, 3-10 s spin periods; 
iHevl fc Hernauistlll998|) . with the absorption reflecting 
either the proton cyclotron line or transit ions between 
bound states of neutral hydro gen (e.g., lHaberl et alJ 
l2003Hvan Kerkwiik et alll2004fl . 

In order to help determine the nature of the ab- 
sorption features, as well as to elucidate what sets 
the isolated neutron stars apart from young rotation- 
powered pulsars, we have started a program to ob- 
tain phase-connected timing solutions, and use these 
to estimate ages and magnetic field strengths. Ear- 
lier, we presented our first resu lts, for RX J0720.4— 3125 
ijKanlan fc van Kerkwiikl|2005l) : here, we consider a sec- 
ond source, RX J1308.6+2127. 

RX J1308.6+2127 (also known as RBS 1223 and 
1RXS J130848.6+212708) was identifie d as a possi- 
ble n earby isolated neutron star by iSchwooe et al.l 
1)1999(1 . The identification was confirmed with the de- 
tection of a 5.16-s X-ray periodicity {H ambarvan et aD 
2002) and a very faint (V w 28 mag) prob- 
able optical counterpart with no radio emission 
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(|Kanlan. Kulkarni. & van Kerkwiikll2002a[K By compar- 
ing the perio ds measured from Chan dra and archival 
ROSAT data, Ha mbarvan et all l)2002|) inferred a spin- 
down rate of P — (0.7 — 2.0) x 10 -11 s s" 1 , implying a 
very strong magnetic field of > 10 14 G. lHaberll {2004), 
however, showed that the 5.16-s periodicity was in fact 
the first harmonic, i mplying a true period of 10.31 s (see 
also lHaberl et al.ll2003l) . and that the spin-down rate in- 
fer red earlier was likely erroneous. This was confirmed 
bv lSchwope et all l|2005|) . who attempted to determine a 
phase-connected timing solution, but failed due to cycle- 
count ambiguities. 

Here, we show that with additional Chandra obser- 
vations specifically obtained for timing purposes, we 
can obtain an unambiguous timing solution. We de- 
scribe our analysis of the Chandra data, as well as 
of archival ROSAT, Chandra, and XMM-Newton data, 
in § |2 and use these to obtain a timing solution in 
§ 01 Since our an alysis closely follows the one u sed for 
RX J0720.4-3125 <|Kaplan fc van KerkwiiU2005|) . we fo- 
cus primarily on those aspects of the analysis that differ. 
We discuss the implications of our result in § fij 

2. OBSERVATIONS 

Our primary data are eight observations taken with 
the Advanced CC D Imaging Spectrometer (ACIS; 
iGarmire et alJl2003f) aboard the Chandra X-ray Obser- 
vatory (CXO). These were designed for timing accu- 
racy, consisting of two sets of four exposures in the 
Continuous-Clocking (CC) mode geometrically spaced 
over a period of about two weeks and separated by about 
half a year. We combined these with data from other 
Chandra observations, as well as from observations with 
XMM-Newton and ROSAT. A log of all observations 3 is 
given in Table 

For the Chandra data, we processed the level- 1 

3 We do not use data from R OS AT observation 70384 8, since 
we, like Hambarvan et al. 12002) and SchwoDe et al. (2005), found 
there were insufficient counts to extract a period or an arrival time. 
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TABLE 1 

Log of Observations and Times of Arrival 



T 

Instrument 


1JJ 


Date 


Exp. 


Counts 


rrn A c 

1 UA 








(ks) 




^IVlJU ) 


HRI 


7040S2 


1998 Oct 01 


4.8 


498 


50824 21 434961' 30 s ) 


ACTS 1 /8 


731 


2000 Jun 24 


9.5 


7395 


51719 5182790(51 


PN/sw 


377-U2 


2001 Dec 31 


18.0 


10633 


52274.2594926(10) 


PN /if 


561-S5 


2003 Tan 01 


27.0 


66219 


52640 4325929("5l 


MOS1 


561-S3 


2003 Jan 01 


29.0 


14469 


52640.4265055(9) 


MOS2 


561-S4 


2003 Jan 01 


29.0 


14925 


52640 4265042(71 


PN/ff 


743-S3 


2003 Dec 30 


30.0 


74587 


53003.4668264(6) 


MOS1 


743-S1 


2003 Dec 30 


32.0 


16086 


53003.4607372(7) 


MOS2 


743-S2 


2003 Dec 30 


32.0 


16523 


53003.4607388(9) 


HRC 


4595 


2004 Mar 30 


90.1 


33823 


53095.3847872(16) 


ACIS CC . . . 


5522 


2005 Feb 14 


16.0 


8953 


53415.6878535(10) 




5523 


2005 Feb 15 


5.7 


3254 


53416.5968822(17) 




5524 


2005 Feb 19 


5.2 


2922 


53420.1717683(14) 




5525 


2005 Mar 10 


5.6 


2923 


53439.0475044(11) 


ACIS CC . . . 


5526 


2005 Jul 09 


15.1 


7684 


53560.2591584(8) 




5527 


2005 Jul 10 


5.1 


2876 


53561.2549594(17) 




5528 


2005 Jul 14 


5.2 


2937 


53565.7697939(18) 




5529 


2005 Jul 29 


5.2 


3048 


53580.7839447(14) 



a HRI: High- Res olution Imag er IZombeck et al.M1995ll aboard ROSAT. PN: XMM-Newton's European Photon Imaging Camera with PN 
detectors IStruder et al,ll2001l) . used in full-frame (i f ) or small window (sw) mode, with thin filter. MOS1/2: European Photon Imaging 
C Cameras wil.U MOS detectors aboard XMM- Newton iTjrrnere^aJ2001|) , used in small-window mode with thin filter. HRC: High- Resolution 
Camera for spectroscopy aboard Chandra (HRC-S: lKraf^et^nH^97ll ! used with the Low-Energy Transmission Crating (LETG). ACIS: 
Chandra's Advanced CCD Imaging Spectrometer, with the S3 chip used either in 1/8-subarray or in Continuous Clocking (CC) mode. 

b Observation identifier (CXO, ROSAT) or revolution number and exposure identifier (XMM-Newton). 

c The TOA is defined as the time of maximum light of the pulsation after the deepest minimum closest to the middle of each observation, 
and is given with 1-<t uncertainties. 



event lists to the level-2 stage following standard pro- 
cedures and the latest calibration set (CALDB ver- 
sion 3.1.0). For the ACIS data, we extracted events 
within 1" of the source with energies between 0.2 and 
2.0 keV, and then applied a clo ck correction of 284.7 /j,s 
Ipavis. Holmes, fc Mversll2003l) . For the HRC-S/LETG 
data, we extracted zeroth-order events from a circle with 
radius of 1.3, and first-order events using the stan- 
dard LETG spectral extraction windows, but limited to 
10 A < A < 65 A. Finally, we used the axbary pro- 
gram to barycenter all of the events (using the X-ray 
position: q. 1? nnn = 1 3 h 08 m 48 / . / 27, 8 J2 ooo = +21°27'06"8; 
Kaplan et al]l2002a|) 

For the XMM-Newton data, we used the standard pro- 
cedures emchain and epchain (XMMSAS version 6.5.0) to 
reprocess the observations. Next, we extracted events 
within 37'.'5 of the source position (using standard qual- 
ity and pattern selections) with energies from 0.12 keV to 
1.2 keV, and used barycen to convert the arrival times to 
the solar-system barycenter. And finally, for the ROSAT 
HRI data, we extracted the events within 12'.'5 of the 
source. We barycentered these using the FT00LS pro- 
grams ABC and BCT, and converted the event times from 
Coordinated Universal Time (UTC) to Barycentric Dy- 
nam i cal T ime (TDB) using the corrections supplied in 
|Co3(|200ll p. 14). 

3. TIMING ANALYSIS 

Our goal is to use times-of-arrival (TOAs) to infer a 
phase-coherent timing solution in which each cycle of 
the source is accounted for. To measure TOAs, we need 
an initial reference period. We det ermined this from th e 
ACIS CC data-sets using a Z\ test ijBuccheri et alJ l983). 
which combines power from the 10.31-s fundamental with 



that from the 5.16-s harmonic (we could not detect sig- 
nificant power in any higher harmonics). We calculated 
Z\ both for the individual sets of four observations (ob- 
servations 5522-5525 and 5526-5529) as well as for the 
eight observations combined. The overall Z\ spectrum 
has a single, well-defined peak, and all possible aliases 
are at considerably lower significances. The measured 
period is P = 10.31252293(9) s (here and below, num- 
bers in parentheses indicate the formal l-cr uncertainties 
in the last digit unless otherwise indicated). 

Using the period derived above, we constructed binned 
light curves (with 16 phase bins) for all of the observa- 
tions. Unlike RX J0720.4-3125, RX J1308.6+2127 has 
non-sinusoidal pulsations and hence determining TOAs 
by fitting a single sinusoid would be inappropriate. In- 
stead we fit for both the first harmonic and fundamental, 
with a possible phase shift between them: 

Ni = A [sin (2tt(M + <f> )) + r 2 sin {^{foU + fa + A<fe))]+C, 

where Ni is the number of counts in bin i, fo — 
0.096969481 Hz is the frequency of the fundamental from 
above, U the time of bin i on the interval [0, P), and the 
parameters are amplitude A, phase (f>o, amplitude ratio 
T2, phase offset A02, and constant offset C. As can be 
seen in Figure ^ this model provides a good fit to even 
the highest-quality lightcurves, which are those derived 
from the long EPIC-PN observations. The inferred phase 
offset and relative amplitude are similar for both obser- 
vations, with A02 w 0.03, and r 2 ~ 2.8. 

Assuming that the pulse profile is not varying, we mea- 
sured arrival times for all of the data using the above 
model, but with Afa and r 2 fixed at 0.03 and 2.8, re- 
spectively, which generally provides a good fit to the 
light curves (see Fig. □ for the EPIC-PN data). Our 
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TABLE 2 

Measured and Derived Timing Parameters for 
RX J1308. 6+2127 From Chandra and XMM-Newton Data 



Pulse Phase 



Fig. 1. — Pulse profiles of the two long EPIC-PN observations, 
showing the different fits to the binned lightcurves. Top: data 
from revolution 561; bottom: data from revolution 743. The solid 
red lines show the fits using Eqn. [Tl with A 02 and free , while 
the dashed blue lines (largely indistinguishable from the red lines) 
show the fits with A</>2 and T2 fixed to the average best-fit values of 
0.03 and 2.8. The zero point in phase has been chosen such that it 
coincides with the maximum following the deepest minimum, which 
we define as the reference position for the TOAs; it is indicated by 
the vertical dashed line at phase 0. 

resulting arrival times are listed in Table |"IJ Here, we 
assigned these times to the time of the maximum of the 
model lightcurve that follows the lowest minimum and is 
nearest to the middle of the observation (see Fig.[U this 
is th e spectrally harder maximum from iSchwone et al.l 
2005). For completeness, we note that the TOAs de- 
rived using the above model are consistent with those 
one would find using a simpler model consisting of only 
the first harmonic (i.e., a sinusoid at P = 5.16 s): the 
mean absolute difference is 0.2cr. 

We determined a timing solu tion from the TOAs us- 
ing an iterative procedure, as in lKanlan fc van Kerkwiikl 
( 2005). We started with the Chandra ACIS CC TOAs 
and found that we were able to fit these without the need 
for a frequency derivative and with no cycle ambiguity, 
consistent with the power spectrum analysis. However, 
between the first ACIS CC TOA and the preceding TOA 
(Chandra HRC) there is a gap of 320 days and we found 
that a constant-frequency model led to a poor fit (x 2 = 
202.7 for 7 degrees of freedom). Including a frequency 
derivative i>, the cycle count over the gap became am- 
biguous by ±1 cycle. The three possibilities — 2683550, 
2683551, and 2683552 cycles — lead to three possible so 
lutions for v: 2.1x10 
-1.1 x 10~ 15 Hz s" 1 
remaining data eliminates the first two of those solutions: 
with the XMM-Newton TOAs (and the remaining ACIS 
point), they have x 2 = 439.1 and 1107.6 with TOA rms 
of 0.041 s and 0.099 s, respectively, while the third has 
X 2 = 18.7 with a TOA rms of 0.010 s (all for 14 degrees 
of freedom) . 

The timing solution derived from the Chandra and 
XMM-Newton TOAs is presented in Table [3 As indi- 
cated by the value of \ 2 , the data are well reproduced by 
simple spin-down; the addition of a cubic (i)) term leads 
to only an insignificant reduction in x 2 , from 18.7 to 18.6, 
and the inferred value of v = —7(15) x 10 -26 Hz s~ 2 is 
consistent with zero. 

Unlike the Chandra and XMM-Newton data, however, 



15 Hz s~\ 1.0xl0~ 15 Hz s" 1 , and 
Fortunately, the addition of the 



Quantity 3 

Dates (MJD) 

t (MJD) 

v (Hz) 

v (Hz s" 1 ) 

TOA rms (s) 

X 2 /DOF 

P(») 

P (s s" 1 ) 

E (erg s _1 ) 

/>'.-. (G) 

Tchar (yr) 



Value 



51720-53581 
53415.687853(2) 
0.0969694896(2) 
-1.053(3) x 10- 15 
0.010 
18.7/14=1.34 

10.31252206(2) 
1.120(3) x 10" 13 
4.0 x 10 30 

3.4 x 10 13 

1.5 x 10 6 



Note. — Uncertainties quoted are twice the formal 1-<t uncer- 
tainties in the fit. 



the ROSAT HRI point is slightly discrepant from the fit 
in Table |2 exceeding it by 0.15(3) cycles. If we include 
it in the fit, the overall solution does not change dras- 
tically (i> becomes -1.050(2) x 10~ 15 Hz s _1 ) and the 
deviation does decrease to 0.11(3) cycles, but the overall 
X 2 increases to 43.5 for 15 degrees of freedom (with a 
TOA rms of 0.0271 s). The deviation could be intrin- 
sic (a glitch, timing noise, precession, changes in pulse 
profile) , but might also be instrumental, perhaps related 
to differences in energy responses between the different 
instruments. Unfortunately, the sampling is too sparse 
to determine a unique solution. As an example, we show 
in Figure |2] a cubic fit that matches the HRI point rea- 
sonably well, with x 2 — 22.5 for 14 degrees of freedom, 
TOA rms of 0.011 s, and v = -2.5(5) x 10~ 25 Hz s~ 2 . 
We stress, however, that, in essence, this solution ap- 
proximates the Chandra and XMM-Newton data with a 
i) = model and then adjusts v to pass through the HRI 
point. Hence, it does not gives much additional insight. 

4. DISCUSSION & CONCLUSIONS 

The spin period and spin-down rate we derive 
for RX J1308. 6+2127 are quite similar to those of 
RX J0720.4— 3125, the only other isolated neutron sta r 
with a timing solution ( Kapla n* fc van Kerkwiikl 12005). 
placing both of them above the pulsar "death-line" in 
a P-P diagram despite t heir lack of radio emission 
l|Kanlan et alJl2002al 12003). Hence, assuming both spin 
down by magnetic dipole radiation, RX J1308. 6+2127 
and RX J0720.4— 3125 have similar inferred magnetic 
field strengths (B dip = 3.4 and 2.4 x 10 13 G, respec- 
tively), characteristic ages (r c har = 1.5 and 1.9 Myr), 
and spin-down luminosities (E — 4.0 and 4.7 x 
10 30 ergs s _1 ). Both sources also have thermal X-ray 
spectra with temperatures of ~ 1 x 10 6 K, but super- 
posed on these are rather different absorption features: 
in RX J1308. 6+2127, the absorption feature is centered 
near or below 300 eV and is very w ide and strong, wit h 
an equivalent width of > 150 eV (|Haberl et alJ 12003): 
in RX J0720.4— 3125, on the other hand, no absorp- 
tion was apparent until 2004, when an absorption fea- 
ture with only a slightly higher energ y but much sm aller 
equivalent width (~ 40 eV) appeared IjHaberl et alJl2004t 
Ide Vries et all200l . If the magnetic fields, line energies, 
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Fig. 2.— Phase residuals for RX J1308. 6+2127. The top panel 
shows the residuals for each TOA compared to a linear (u = 0) 
model. The solid curve gives the best-fit quadratic (v 7^ 0, v = 0) 
ephemeris for all the Chandra and XMM-Newton data (Tab. 0. 
The vertical dashed line indicates the reference time to- The bot- 
tom panel shows the residuals relative to the quadratic model; the 
ROSAT HRI point is off the scale at 0.14 cycles, and so is indicated 
with an arrow in this panel. We also show a best-fit cubic model 
fitted to all of the data including the ROSAT point (i> 7^ 0; dotted 
line). 



and blackbody temperatures are all similar, what would 
account for the significant difference in line strength? 

The interpretation of the absorption li nes is still a 
matte r of debate. For RX J1308.6+2127. lHaberl et all 
(2003) suggested proton cyclotron absorption in a field 
of (2-6) xlO 13 G, which could match both the energy and 
strength of the observed feature (for a gravitational red- 
shift of z = 0.3, the observed line energy would be at 
100-300 eV). As the feature in RX J 0720.4-3125 was so 
much weaker but at a similar energy. Ivan Kerkwiik et al.l 
l)2004[) argued that for that source the proton cyclotron 
line was below the observed band and that the absorption 
was due to the — > 2 transition between tightly-bound 
states of neutral hydrogen in a ~2 x 10 13 G field. 

Qualitatively, the field strengths inferred from our tim- 
ing measurements agree with the above expectations: 
the field of RX J1308.6+2127 is stronger than that of 
RX J0720.4— 3125, while a weaker field would be ex- 
pected if the absorption in both sources were due to 
the same mechanism. Quantitatively, a larger difference 
between the two sources was expected. However, this 
discrepancy may simply reflect orientation and/or sub- 
structure in the magnetic field — higher-order multipoles 
near the surface would affect the emission properties but 
not the spin-down rate. Such effects may also be respon- 
sible for the difference in pulse profile: double-peaked for 
RX J1308.6+2127, and sinusoidal for RX J0720.4-3125. 
Phase-resolved X-ray spectroscopy coupled with obser- 
vations of more sources are probably the best ways to 
disentangle these effects, and efforts are underway (e.g., 
lHaberl et ai]l2004USchwope et al.ll2005[) . 

Another open issue is that the characteristic ages 
of RX J0720.4-3125 and now RX J1308. 6+2127 are 
three to four times larger than the values of ~ 
0.5 Myr one would expect for simple cooling models 
llHevl fc HernauistllQMlKaplan et al.l2002bHZane et all 
120021) and, in the case of RX J0720.4-3125, tracing the 
object back along its trajectory to a likely birt h location 
ijMotch. Zavlin. fc Haberll 12001 lKapIar] |2004). Indeed, 



this age discrepancy may extend to the other isolated 
neutron stars discovered by ROSAT: all have similar 
temperatures and thus likely similar cooling ages, and 
most also have similar periods (3-10 s) and, based on 
the similar energies at which they show X-ray absorption 
features (0.3-0.7 keV), similar magnetic field strengths, 
implying similar characteristic ages. 

We first consider whether the discrepancy could 
result from the characteristic age being an over- 
estimate. In general, for a spin-down torque cx 
v n , the pulsar's spin-down age is given by t s d = 

P/{n - 1)P 1 - (Po/P)" -1 , where P is the initial 

spin period and n = vvjv 1 is the "braking index," 
equal to 3 un der the assumption of magn etic dipole ra- 
diation (e.g., iManchester fc Tavlorl 119771 p. 111). For 
Pq <C P and n — 3, one recovers the characteristic age 
T c har = P/2P, but if Pq is not much smaller than the cur- 
rent period P, the characteristic age is an overestimate. 
While we cannot exclude the required birth periods of 
7-8 s, there is as yet no concrete indication that neutron 
stars are born with periods longer than 1 00 ms (e.g., 
IKasni fc Helfanll27fflalGotthelf et alJl200l . 

For RX J0720.4-3125, this led us to discuss the pos- 
sibility that the neutron star formed in a binary sys- 
tem and accreted matter before being ejected in a sec- 
ond supernova; at ejection, the neutron star would 
still be hot from the accretio n, but spinning slowly 
ijKaplan fc van Kerkwiik! 12005(1 . A variation on this 
model would involve accretion from a residual debris disk 
such as that recently disc overed around the anomalous X- 
ray p ulsar 4U 0142+61 (Wang, Ch akrabartv. fc Kaplan] 
2005), but this assumes that the accretion disk persists 
for a sufficient time and can affect the spin-down of the 
neutron stars, both of which are far from clear. In ei- 
ther case, if accretion played a role, it might explain why 
the spectral properties of the isolated neutron stars are 
rather different from those of the radio pulsar population. 

Unfortunately, the above solution for the age discrep- 
ancy is not unique. First, it is possible that the spin- 
down was not due to a constant dipole but that the 
magnetic field decayed (effectively, this implies n > 3). 
Second, the cooling age could be incorrect due to non- 
standard cooling. For instance, the energy released by 
mag netic field decay migh t keep a neutron star hot- 
ter iJHevl fc Kulkarnil [T998: note, however, that these 
authors do not expect significant effects for the field 
strengths we infer) and longer cooling times are also ex- 
pected for a light neutron star ijYakovlev et alJl2004j) . A 
prolonged cooling timescale seems somewhat unlikely in 
light of the kinematic age estimates, though it should 
be kept in mind that these may also not be unique: for 
RX J 0720.4— 3125, several different birth places are pos- 
sible <|Motch et alJl25o3 lKapladl2004|) . 

Finally, when fitted with just a simple spin-down 
model, RX J0720.4— 3125 showed significant timing 
residuals of 0.31 s (root-mean-square), much larger 
than the uncertainties on the TOAs. In contrast, for 
RX J1308. 6+2127 the root-mean-square residuals are 
only 0.010 s (excluding the HRI point), consistent with 
measurement errors. Since the sources were observed 
with the same instruments, and have similar spectral 
shapes, the differences in timing behaviors are likely not 
instrumental. Instead, they may well reflect the fact that 
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both the spectrum and pulse profile of RX J0720.4— 3125 We acknowledge support through Chandra grant G05- 

are v arying with time (|de Vries et all 120041: iVink et all 6050 A. 
|2004|) . while RX J1308.6+2127 appears to be stable. 
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